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Microcanonical analysis is a powerful method for studying phase transitions of finite-size systems. 
This method has been used so far only for studying phase transitions of equilibrium systems, which 
can be described by microcanonical entropy. I show that it is possible to perform microcanonical 
analysis of a nonequilibrium phase transition, by generalizing the concept of microcanonical entropy. 
One-dimensional asymmetric diffusion process is studied as an example where such a generalized 
entropy can be explicitly found, and the microcanonical method is used to analyze a nonequilibrium 
phase transition of a finite-size system. 
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Microcanonical analysis is a powerful method for 
studying phase transitions of finite-size systems 0,0. 
In this approach, the form of the microcanonical entropy 
is examined to see whether there is a convex region. Ex¬ 
istence of such a region signals the onset of an inhomo¬ 
geneity, and the system is considered to undergo a first- 
order phase transition in this region. The phase tran¬ 
sition in microcanonical analysis is a well-defined con¬ 
cept even for a finite-size system, which is in contrast 
to most of the traditional canonical ensemble approach 
where phase transition is defined only for infinite size 
systems, defined in terms of singular behavior of physi¬ 
cal quantities in this limit. The microcanonical analysis 
has been applied for studying phase transitions of vari¬ 
ous finite-size systems such as spin models , atomic 
clusters and nucl ei 0 , 0 , polymers [m, peptides [l6| . 
and proteins 

Since only equilibrium systems can be described in 
terms of microcanonical entropy, microcanonical method 
has been used exclusively for analyzing equilibrium phase 
transitions so far. In this Letter, I show that it is possi¬ 
ble to apply this method for analyzing a nonequilibrium 
phase transition 22 h 28 I|. by a proper generalization of 
the concept of microcanonical entropy. 

Let us first briefly review the connection between the 
convex region of microcanonical entropy and the phase 
transition 0,0 . We consider a finite closed system with 
a conserved quantity, say energy E, and denote the num¬ 
ber of corresponding microstates as UIl{E), where the 
subscript denotes the dependence on the system size L. 
The microcanonical entropy is then defined as 


SL{E)=\nnL{E) 


( 1 ) 


where we use the unit with ks = I- Now suppose we 
construct a larger system by assembling two identical 
subsystems of energy E and size L. We let the two 
subsystems make a thermal contact, but let the cou¬ 
pling between the two systems be weak enough so that 
the total energy is Etot = 2E = Ea + Eb where Ea 
and Eb are the energy values of the two subsystems. 
We then examine the qualitative feature of the probabil¬ 


ity distribution of the energy values of the subsystems, 
P{Ea,Eb) (xe:>ip{SLiEA) + SLiEB)). li Sl^E) is a con¬ 
cave function, then P{E,E) > P{E-,Ea) for any E± 
with E- < E < E+, so the homogeneous distribution of 
energy among the subsystems is preferred. On the other 
hand, if there is a convex region in Sl{E) so that one 
can find values Ex , E 2 and 0 < p < I satisfying 

E = pEi + {l—p)E2, 

SUE) < pSUEi) + il-p)SUE2), ( 2 ) 

then there are values E± such that P{E, E) < 
P{E-, EU so that an inhomogeneous distributions is fa¬ 
vored, and we say that the system is in the region of 
the first-order phase transition. The argument can eas¬ 
ily be generalized to the case of subsystems of different 
sizes, more than two subsystems, and multiple conserved 
quantities 00. 

We note that the only relevant property of UIl{E) ex¬ 
ploited in the argument is that when a conserved quantity 
Q = Qa -}- Qb of the total system is distributed over two 
subsystems A and B, the probability distribution of Qa 
and Qb is proportional to the product of fls: 

^’(Qa, Qb) oc ni,(QA)fJL'(QB) (3) 

where L and L' denote the sizes of the subsystems A and 
B. We can also impose a certain boundary condition 
at the interface between the subsystems, in which case 
^’(QajQb) becomes a conditional probability. There¬ 
fore, it is clear that even for a nonequilibrium system, if 
a probability of the distribution of a conserved quantity 
Q among subsystems under appropriate boundary con¬ 
dition can be expressed in the form Eq.ljS]), then we can 
consider nL(Q) as the generalized density, and their log 
as the generalized entropy, which can then be used as the 
target of the microcanonical analysis. This is the main 
claim of this Letter. 

As an example of a nonequilibrium model on which 
microcanonical analysis can be performed, we consider a 
diffusion model where particles of two types, labelled as I 
and 2, move asymmetrically on a periodic lattice of length 
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L 25, [ 2 ^ . Treating the vacancy as a particle with 

label 0, the transition rates gap for the particle exchange 
of the type (a, (3) —>■ (/3, a) at neighboring sites are given 


510 = 502 = 1, 512 = 5, 521 = 1 (4) 


with all other components of gs being zero. We note that 
the numbers of both types of particles are conserved sep¬ 
arately, which we will denote as rii and 712 - The matrix 
representation of the stationary state for this process has 


already been found, and is given as [2J, [25 


PstiPi, • • • Pl) — trG/3j ■ • ■ 


E 


G^,(5) 


where Pk denotes the particle type at the fc-the site, and 
the components of the three infinite-dimensional matrices 
Gp (/3 = 0,1, 2) are given as 

(Go)y = SiiSij 
(Gi)jj — OjiSij 

(G2)ij = QiSij + SjSi-ij. (6) 

where 

1-H - 2gi-fe 
5 - 1 

Now let us suppose that there are vacancies at sites a 
and b. The whole periodic lattice can be divided into two 
regions bounded by these two sites, and we would like to 
obtain conditional probability for particles in these two 
regions being ua = and 

Obviously, from Eq.([5l), we see that it is proportional to 


a-k = 

Sktk = 


P{nA,nB) oc 

E ■ ■ 


,_iGoG-y^^i • • • 


7i)"^7L 



xGoG^,,_^j • • 



)^(E<57,.2,4^^) 




j&A 

x^(E 


nf'^) 

1 (8) 

k^B 

1GB 




where 5{a,b) = Sa,b denotes Kronecker delta function 
that vanishes when the indices are not equal. Note that 


tr(G^i ■ • • G-y^_iGoG-y^^i ■ • • G^^_jGoG-yj^j • • ■ G-^^) 
= tr(GoG^^^j • • • G^j_j GoG^,,_^j • • • G^^_ J 

= [*^7a + l ■ ■ ■ J • • • G-y^_J 

= tr(GoG-y^^j • • • G^,,_ Jtr(GoG^,,^^ • • • G^^_ J. (9) 

Therefore, the conditional probability for the steady state 
is expressed in the form Eq.(|31), where the generalized 


density for a system of size L is now defined as 


Oi,(n) 


L-l 


L-1 


L-l 


= E tr(Go n ^'r'=)^(E'^'^•4 >^i)^(E 


71O-07L 


/c=l 




= E 


L-l 

.fe=l 


7fc 


L-l L-l 

^(E '^7i.l>^l)^(E '57j,2,«2)(10) 
11 i=l 1 = 1 


where the system size L includes one vacancy. We an¬ 
alyze the phase transition of the current model by per¬ 
forming the microcanonical analysis on the generalized 
entropy 5'L(n) = logni:,(n). It is expressed in terms of 
(L/2) X (T/2) submatrices of Gp, which can be computed 
exactly for given values of q and L 24, 25, 2^. 

By performing analytic computations, Monte Carlo 
simulation, mean field calculations 2^, and partition 
function zero analysis [^ . it has been argued that this 
system undergoes a nonequilibrium phase transition in 
the limit of L —> 00 . There is a gc > 1 such that the 
system remains homogeneous for g > gc, but inhomo¬ 
geneities of particle densities appear for certain range of 
particle numbers when q < q^. In fact, the latter can be 
considered as a region of first order transition between 
the fluid and the condensed phases, as will be elaborated 
below. 

From the viewpoint of microcanonical analysis, the cri¬ 
terion for a first-order transition is the existence of a 
nonconcave region in the microcanonical entropy, a set 
of points where one can find a direction with positive 
second derivative Hi- For the current model where the 
conserved quantity n is discrete, I examined discretized 
second derivatives 


AaAbSL = SL{ni + a,n2 + b) + SL{ni - a,n2 - b) 

- 25 ^( 711 , 712 ), (II) 

along the horizontal ((a, &) = (1,0)), vertical ((a, &) = 
(0,1)), and two diagonal ((a, 6) = (1,±1)) directions. A 
point in the interior is nonconcave if any one of these four 
quantities has a positive value. 

The generalized entropy function 5 l(71i, 712 ) are shown 
in the left panels of the figures [I] and [5] for L = 10 and 
L = 100 respectively, for various values of g. Note that 
the entropy has the symmetry with respect to the line 
Til = TI 2 due to the invariance under the simultaneous 
application of particle type exchange 1 o 2 and the par¬ 
ity inversion k ^ —k. We see that for small enough 
values of g, a nonconcave region appears in the gener¬ 
alized entropy, enclosed by dashed lines in figure |T] and 
denoted as gray regions in figure [H As g increases, the 
nonconcave region shrinks, and eventually disappears for 
large enough values of g. We find that nonconcave re¬ 
gion always includes a part of the line rii = 712 . The 
second derivative at such a point is also largest along the 
(1,1) direction, which tells us that for sufficiently small g. 


















3 



FIG. 1: The generalized entropy function SL{ri\,n 2 ) for L = 
10 is displayed with distinct symbols for points belonging to 
different ranges of function values, for (a) q = 0.5 , (c) q = 0.7, 
and (e) q = 2.0. The nonconcave region is enclosed by dashed 
lines. The cross section along the diagonal line ni = n 2 in 
figures (a), (c), and (e), are displayed in figures (b), (d), and 
(f). The concave envelopes are drawn in figures (b) and (d) 
with dashed lines as visual guides. 


FIG. 2: The contours of the generalized entropy function 
SL{ri\,n 2 ) for L = 100 are drawn for (a) q = 0.8 , (c) q — 1.1, 
and (e) q = 2.0, at intervals of 200, 50, 20, respectively. The 
nonconcave region is colored gray. The cross section along the 
diagonal line ni = n 2 in figures (a), (c), and (e), are displayed 
in figures (b), (d), and (f). The concave envelopes are drawn 
in figures (b) and (d) with dashed lines as visual guides. 


when the system is divided into subsystems with respect 
to a pair of vacancies, it is most probable that there is a 
inhomogeneity for the total particle numbers, but there 
are the same numbers of two species at both sides. Note 
that this is an exact statement for a finite value of L, in 
contrast to the results of previous works where limit of 
L —oo was considered [^, . 

The cross sections of SL{ni,n 2 ) along the line ni = 
712 = n, SL(n, n), are also displayed in the right panels of 
figures [Hand [21 where the concave envelopes are denoted 
by dashed lines whenever they exist. These correspond 
to the regions of the first-order transition, whose upper 
and lower boundaries p± in the space of particle density 


p = n/L {0 < p < 0.5) are drawn as functions of q to 
produce a phase diagram in the figure [31 for L = 10 and 
L — 100. The mean field result q{p) = (1 -b 6/9)/(l -b 2p) 
in the limit of L —^ oo is shown in the figure [31 with a 
dashed line for comparison , where q{p) is the inverse 
function of p±{q). 

The high-density side of the phase boundary, p > p^, 
corresponds to the condensed phase. Note that there is 
a qi{L) such that p- = 0 for q < qi{L), in which case 
the low-density phase p = p- = 0 is just the vacuum 
without any particles present. For q > qi{L), the low- 
density phase p < p- is the fluid phase. As q increases, 
p± approaches toward each other and eventually merges 
at the critical point q = qdL), after which the system 
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FIG. 3: The phase boundary p±{q) for L = 10 and L = 100. 
The mean field result in the limit of L —>■ oo is shown with 
dashed line for comparison. 

is in a homogeneous phase. The values of qi and Qc for 
L = 10 and L = 100 are indicated by arrows in Figure 
[3l The mean field prediction for these parameters are 
gi(oo) = 1 and qc{oo) = 2, as can be easily read off from 
the analytic expression for q{p). 

The regions q < qi, qi < q < qc, and q > qr have been 
called pure, mixed, and disordered phases [2^. However, 
microcanonical analysis shows that in p space, each of the 
regions q < qc and qi < q < qc is divided into vacuum 
(or fluid) phase {p < p_), condensed phase (p > p+), 
and the phase coexistence region (p_ < p < p^). The 
situation is analogous to the two-dimensional Ising model 
with the conserved magnetization M and the tempera¬ 
ture T. When one simply considers the T dependence, 
then there is a critical temperature Tc such that the sys¬ 
tem is in a disordered phase for T >Tc and ordered phase 
for T < Tc- However, by examining the M dependent be¬ 
havior of the system, one realizes that the ordered phase 
in fact gets divided into up-spin phase, down-spin phase, 
and the region of the first-order transition between up 
and down phases. 

I also plot qc(,L) and qi{L) as functions of L in figure 
m Both qc{L) and qi{L) approach their mean field values 
qc(po) = 2 and qi(oo) = 1. 

The current work also clarifies the physical meaning 
of a previous work based on the partition function ze¬ 
ros (PFZs) [ 2 ^. There, a partition function of the form 
Sni rt 2 was constructed where x was 

called the chemical potential. Then the PFZs in the 
complex plane of x was analyzed to claim that there is 
a first-order transition as L —>■ 00 , for sufficiently small 
values of q. It is obvious that Hi(n) was used implic¬ 
itly as the generalized density of states, but it was not 
explained why Hi(n) should have such a special status. 
Also, the physical meaning of the chemical potential was 



FIG. 4: The critical values qc (solid line) and qi (dashed line) 
as functions of the system size L. The first-order phase tran¬ 
sition exists for q < qc- The low-density phase is a vacuum 
phase for q < qi and fluid phase for qi < q < qc 


unclear, because Hi(n) was regarded as describing the 
the particles on a periodic lattice of size L, which is an 
isolated system. The current work not only justifies the 
use of Hi(n) as a generalized density, via the factoriza¬ 
tion Eq.Q, but also shows that fli(n) in PFZs approach 
describes a subsystem of size L — 1 bounded by pair of va¬ 
cancies, rather than the whole system. Then the chemical 
potential x can be considered as a parameter describing 
the rest of the system whose size is much larger than L, 
acting as an infinite-size particle reservoir. The micro- 
canonical analysis is more general since the phase transi¬ 
tion is well defined for a system with a finite size. In fact, 
the notion of finite-size nonequilibrium phase transition 
itself is introduced for the first time in the current work 
via microcanonical analysis, which would be a subject of 
much interest for future study. 
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